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Abstract. We present a summary of gamma-ray millisecond pulsar (MSP) observations with the 
Fermi Large Area Telescope. The radio and gamma-ray light curves of these MSPs have been 
modeled in the framework of the retarded vacuum dipole magnetic field. Likelihood fitting of 
the radio and gamma-ray light curves with geometric emission models allows us to give model- 
dependent confidence contours for the viewing geometry in these systems which are complementary 
to those from polarization measurements. 

Keywords: gamma rays: observations, pulsars: millisecond 
PACS: 95 



1. INTRODUCTION 

Soon after the first millisecond pulsar (MSP) was discovered [6] there was much specu- 
lation regarding high-energy (HE, >0.1 GeV) MSP emission (e.g., [22]). Most models 
assume that pulsar HE emission is the result of curvature radiation from electrons and/or 
positrons. In these models the particles are accelerated either near the stellar surface 
(above the magnetic polar cap, e.g. [10]) or in the outer magnetosphere (up to the light 
cylinder) in narrow accelerating gaps with either a two-pole caustic (TPC, e.g. [11]) 
or outer-gap (OG, e.g. [8]) geometry. The derived magnetic fields of MSPs, assuming 
dipole spin down, lie below the curvature radiation pair death line [15] on the P-P dia- 
gram and thus narrow accelerating gaps were not expected. This led to the development 
of the pair-starved polar cap (PSPC) model in which the full open field line volume 
above the polar cap is available to accelerate particles [19]. 

Upper limits (3a) were calculated for 19 MSPs in EGRET data [12] and a 4.9c pulsed 
detection using EGRET data was reported for PSR J0218+4232 [17]. Recently, 4.2cr 
pulsations were detected from PSR B 1 82 1 —24 in the globular cluster M28 using AGILE 
data [20]. It was not clear if these singular cases were unique or if HE emission was 
common in MSPs until the launch of the Fermi Gamma-ray Space Telescope {Fermi). 



2. MILLISECOND PULSARS AND FERMI 



The main instrument aboard Fermi is the Large Area Telescope (LAT) [5], a pair- 
conversion telescope sensitive to gamma rays with energies from 0.02 to >300 GeV. 
Using ~8 months of data, significant HE pulsations from 8 MSPs were detected with the 
LAT [1]. Some MSP gamma-ray light curves showed sharp peaks reminiscent of what 
is seen in younger gamma-ray pulsars [2]. This implies that the HE emission processes 
in MSPs must be the same as those in younger gamma-ray pulsars. In particular, MSPs 
are able to form narrow accelerating gaps despite being below the theoretical curvature 
radiation pair death line [15]. This was supported by modeling of the gamma-ray and 
radio light curves for these 8 MSPs [23] using TPC, OG, and PSPC models for the 
gamma-rays and a hollow-cone beam for the radio [21]. It was found that 6 were well fit 
with outer-magnetospheric models while only 2 required the PSPC model. 

Pulsed gamma rays have also been detected from PSR J0034— 0534 in which the 
gamma-ray and radio profiles are aligned [3], a phenomenon previously seen only 
in the Crab pulsar. To match the observations, the radio emisson was modeled as 
significantly extended in altitude, near the light cylinder, and contained within the 
gamma-ray emisson region. 

In addition to the pulsed detection of individual MSPs, gamma-ray point sources have 
been found consistent with 8 globular clusters [4] known or suspected of harboring many 
MSPs, which also display the characteristic pulsar spectrum, implying that the observed 
emission is from a combination of many MSPs. Radio searches of unassociated Fermi 
LAT sources have revealed 23 new MSPs to date, some of which have now been seen to 
pulse in gamma-rays as well (e.g., [14]). 

3. SIMULATIONS AND LIKELIHOOD FITTING 

We have simulated radio and gamma-ray light curves following the procedures of [23], 
using the retarded vacuum dipole magnetic field geometry, with resolutions of 1° in 
inclination angle (a), the same in viewing angle (£), and 0.05 in gap width (normalized 
to the polar cap opening angle). Electrons are followed along the magnetic field lines out 
to a radial distance of 1.2 Rlc but not beyond a cylindrical distance of 0.95 Rlc> where 
Rlc = cP/(2tt). We have also included the Lorentz transformation of the magnetic field 
from the inertial observer's frame to the co-rotating frame which has been advocated as 
necessary for self-consistency [7]. To fit the light curves we have developed a Markov 
chain Monte Carlo (MCMC) maximum likelihood procedure which uses small-world 
chain steps [13] and simulated annealing [18]. An MCMC involves taking random steps 
in parameter space and accepting a step based on the likelihood ratio with respect to 
the previous step [16]. It is necessary to balance the precision of the observed radio 
profiles against the simplistic cone-beam geometry we use and the best choice of radio 
uncertainty is under investigation. For the example presented in Section 4 we use an 
uncertainty equal to the average relative gamma-ray error, in the on-peak region, times 
the maximum radio value. 



4. PRELIMINARY RESULTS 



Our MCMC algorithm has been applied to PSR J2017+0603, a 2.896 ms pulsar which 
is one of three MSPs discovered to date with the Nancay radio telescope in searches of 
Fermi LAT unassociated sources [9, 14]. We find a = 16° and £ = 68° with an infinitely 
thin gap for the TPC model and a = 36° and £ = 74° with a gap width of 0.05 for the OG 
model, see Fig. 1. For the TPC model, an infinitely thin gap is unphysical which implies 
that the best-fit value is less than our width resolution of 0.05. The observed gamma-ray 
light curve is reproduced well by both models with the OG model slightly preferred by 
the likelihood. Both models reproduce the correct radio-to-gamma phase lag but neither 
is successful at generating all of the observed radio features. The complexity of the radio 
profile suggests that the emission may be from a region extended in altitude or that the 
radio beam shape is much more complex than a single- altitude cone. 




FIGURE 1. Data and best-fit light curves for PSR J2017+0603 with gamma-ray (a) and radio (b). TPC 
models are the large dashed lines while OG models are the small dashed lines. 

By examining the distribution of a-£ pairs in the output chains it is possible to obtain 
confidence contours in viewing geometry marginalized over the other parameters (Fig. 
2 for the TPC model). The resulting contours are not constraining but will improve with 
more gamma-ray data. Examining the marginal distributions of the other fit parameters 
and a reveals very asymmetric probabilities which leads to the offset of the a-£ peak 
and the best-fit geometry. There are currently no geometric constraints from radio 
polarization measurements of PSR J2017+0603 for comparison. 
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FIGURE 2. Marginalized confidence contours for the TPC fit of PSR J20 17+0603: within the darkest 
grey area is 39% confidence, next is 68%, and lightest grey area is 95%. The best-fit geometry is indicated 
by the dashed lines. 




5. CONCLUSIONS 



MSPs have been established as a class of HE emitters and observations with the LAT 
have led to a reassessment of MSP HE emission models. We have used geometric models 
to simulate the radio and gamma-ray light curves of PSR J20 17+0603 and produced 
confidence contours in viewing geometry. We plan to apply this technique to all LAT 
detected MSPs. 
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